5th Workshop "Group Analysis of Differential Equations & Integrable Systems" 2010, [T HT31 



Point Symmetry Group 
of the Barotropic Vorticity Equation 

Alexander BIHLO t and Roman O. POPOVYCH^ 

' Fakultat fur Mathematik, Universitat Wien, 

Nordbergstrafie 15, A-1090 Wien, Austria 

E-mail: alexander.bihlo@univie.ac. at, rop@imath.kiev.ua 
j — i 

* Institute of Mathematics of NAS of Ukraine, 
3 Tereshchenkivska Str., Kyiv-4, Ukraine 

The complete point symmetry group of the barotropic vorticity equation on the 
00 . /3-plane is computed using the direct method supplemented with two different 

techniques. The first technique is based on the preservation of any megaideal 
of the maximal Lie invariance algebra of a differential equation by the push- 
forwards of point symmetries of the same equation. The second technique 
involves a priori knowledge on normalization properties of a class of differ- 
'■^ | ential equations containing the equation under consideration. Both of these 

C$ i techniques are briefly outlined. 



Oh 



On 



1 Introduction 

CO 

It is well known that it is much easier to determine the continuous part of the 
complete point symmetry group of a differential equation than the entire group 
including discrete symmetries. The computation of continuous (Lie) symmetries 
is possible using infinitesimal techniques, which amounts to solving an overdeter- 
mined system of linear partial differential equations (referred to as determining 
equations) for coefficients of vector fields generating one-parameter Lie symmetry 
groups. Owing to the algorithmic nature of this problem, the automatic computa- 
^ ■ tion of Lie symmetries is already implemented in a number of symbolic calculation 

packages, see, e.g., papers 0E1G9] for detail description of certain packages and 
reviews [6|[TTJ] . 

The relative simplicity of finding Lie symmetries of differential equations is 
also a primary reason why the overwhelming part of research on symmetries is 
devoted to symmetries of this kind. See, e.g., the textbooks [Hl5lfT9ti2T] for general 
theory and numerous examples and additionally the works [TH3ll8l[T8] for several 
applications of Lie methods in hydrodynamics and meteorology. 

As continuous symmetries, also discrete symmetries are of practical relevance 
in a number of fields such as dynamical system theory, quantum mechanics, crys- 
tallography and solid state physics. They can also be helpful in some issues related 
to Lie symmetries, e.g. allowing for a simplification of optimal lists of inequivalent 
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subalgebras, and due to enabling the construction of new solutions of differential 
equations from known ones. It is not possible, in general, to determine the whole 
point symmetry group in terms of finite transformations by usage of infinitesi- 
mal techniques. On the other hand, the direct computation of point symmetries 
based on their definition boils down to solving a cumbersome nonlinear system 
of determining equations, which is difficult to be integrated. Similar determining 
equations also arise under calculations of equivalence groups and sets of admissi- 
ble transformations of classes of differential equations by means of employing the 
direct method. In order to simplify the derivation of the determining equations, 
different special techniques have been developed involving, in particular, the im- 
plicit representation of unknown functions, the combined splitting with respect 
to old and new variables and the inverse expression of old derivative via new 

ones P3H23I2B]. 

There exist two particular techniques that can be applied for a priori sim- 
plification of calculations concerning the point symmetry groups of differential 
equations. 

The first technique was presented in [11] for equations whose maximal Lie 
invariance algebras are finite dimensional. It is based on the fact that the push- 
forwards of point symmetries of a given system of differential equations to vector 
fields on the space of dependent and independent variables are automorphisms 
of the maximal Lie invariance algebra of the same system. This condition yields 
restrictions for those point transformations that can qualify as symmetries of the 
system of differential equations under consideration. We will adopt this technique 
to the infinite dimensional case using the notion of megaideals of Lie algebras, 
which are the most invariant algebraic structures. 

The second technique involves available information on the set of admissible 
transformations of a class of differential equations [25] , which contains the inves- 
tigated equation. 

In the present paper, we will demonstrate both of these techniques by com- 
puting the complete point symmetry group of the barotropic vorticity equation 
on the /3-plane. This is one of the most classical models which are used in geo- 
physical fluid dynamics. The techniques to be employed are briefly described in 
Section [21 The actual computations using the method based on the corresponding 
Lie invariance algebra and that involving a priori knowledge on admissible trans- 
formations of a class of generalized vorticity equations are presented in Section [3] 
and[H respectively. A short summary concludes the paper. 

2 Techniques of calculation 

of complete point symmetry groups 

Both the techniques described in this section should be considered merely as tools 
for deriving preliminary restrictions on point symmetries. In either case, calcula- 
tions must finally be carried out within the framework of the direct approach. 
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2.1 Using megaideals of Lie invar iance algebra 

The most refined version of the technique involving Lie symmetries in the calcu- 
lations of complete point symmetry groups was applied in [11]. It is outlined as 
follows: Given a system of differential equations C whose maximal Lie invariance 
algebra g is n-dimensional with a basis {e\, . . . , e n }, n < oo, one has to compute 
the entire automorphism group of g, Aut(g). Supposing that T is a transfor- 
mation from the complete point symmetry group G of £, one has the condition 
%ej = Yl7=l e i a ij f° r J = l)---i n ) where % denotes the push- forward of vector 
fields induced by T and (ajj) is the matrix of an automorphism of g in the chosen 
basis. This condition implies constraints on the transformation T which are then 
taken into account in further calculations with the direct method. 

The method we propose here is different to those described in the previous 
paragraph. In fact, it uses only the minimal information on the automorphism 
group Aut(g) in the form of a set of megaideals of g. Due to this, it is applicable 
also in the case when the maximal Lie invariance algebra is infinite dimensional. 
The notion of megaideals was introduced in [24J . 

Definition 1. A megaideal i is a vector subspace of g that is invariant under any 
transformation from the automorphism group Aut(g) of g. 

That is, we have Ti = i for a megaideal i of g, whenever T is a transformation 
from Aut(g). Any megaideal of g is an ideal and characteristic ideal of g. Both 
the improper subalgebras of g (the zero subspace and g itself) are megaideals of g. 
The following assertions are obvious. 

Proposition 1. If\\ and\2 are megaideals of q then so are ii+i2, iiHi2 and [ii , 12] j 
i.e., sums, intersections and Lie products of megaideals are again megaideals. 

Proposition 2. If 12 is a megaideal of \\ and \\ is a megaideal of g then \2 is a 
megaideal of g, i.e., megaideals of megaideals are also megaideals. 

Corollary 1. All elements of the derived, upper and lower central series of a Lie 
algebra are its megaideals. In particular, the center and the derivative of a Lie 
algebra are its megaideals. 

Corollary 2. The radical r and nil-radical n (i.e., the maximal solvable and nilpo- 
tent ideals, respectively) of g as well as different Lie products, sums and intersec- 
tions involving q, x andn ([q,x], [t,t], [fl,n], [r, n], [n, n], etc. ) are megaideals of g. 

Suppose that g is finite dimensional and possesses a megaideal i which, without 
loss of generality, can be assumed to be spanned by the first k basis elements, 
i = (ei, . . . , e&). Then the matrix ((%•) of any automorphism of g has block 
structure, namely, ajj = for i > k. In other words, in the finite dimensional 
case we take into account only the block structure of automorphism matrices. 
This is reasonable as the entire automorphism group Aut(g) (which should be 
computed within the method from |11| ) may be much wider than the group of 



4 



A. Bihlo and R.O. Popovych 



automorphisms of q induced by elements of the point symmetry group G of C. 
Moreover, it seems difficult to find the entire group Aut(g) if the algebra q is 
infinite dimensional. At the same time, in view of the above assertions it is easy 
to determine a set of megaideals for any Lie algebra. 

2.2 Direct method and admissible transformations 

The initial point of the second technique is to consider a given pth order system C° 
of I differential equations for m unknown functions u = (u , . . . , u m ) of n indepen- 
dent variables element of a class C\s of similar systems £$: 
L(x, 9(x, = parameterized by a tuple of pth. order differential func- 
tions (arbitrary elements) = (9 1 (x,u^ p j), . . . ,9 k {x,u^)). Here denotes the 
set of all the derivatives of u with respect to x of order not greater than p, includ- 
ing u as the derivatives of order zero. The class C\s is determined by two objects: 
the tuple L = (L 1 , . . . ,L l ) of I fixed functions depending on x, and 9 and 9 
running through the set S. Within the framework of symmetry analysis of differ- 
ential equations, the set S is defined as the set of solutions of an auxiliary system 
consisting of a subsystem S{x, u^,9^(x, u^)) = of differential equations with 
respect to 9 and a non- vanish condition T,(x,u^,9^(x,u^)) ^ with another 
differential function S of 9. In the auxiliary system, x and it( p ) play the role of 
independent variables and 9^ stands for the set of all the partial derivatives of 
9 of order not greater than q with respect to the variables x and U( p y In view of 
the purpose of our consideration we should have that C° = Ce Q for some 9q £ S. 

Following [25], for 9,9 E S we denote by T(9,9) the set of point transforma- 
tions which map the system Cg to the system Cg. The maximal point symmetry 
group Gq of the system Cg coincides with T(9,9). 

Definition 2. T(£| 5 ) = {(0,0, ¥>) \ 9,9 £ S,<p £ T(9,9)} is called the set of 
admissible transformations in £\$. 

Sets of admissible transformations were first systematically described by King- 
ston and Sophocleous for a class of generalized Burgers equations [H] and Winter- 
nitz and Gazeau for a class of variable coefficient Korteweg-de Vries equations [30J , 
in terms of form-preserving |14H16| and allowed |30| transformations, respectively. 
The notion of admissible transformations can be considered as a formalization of 
their approaches. 

Any point symmetry transformation of an equation Cg from the class £\s gen- 
erates an admissible transformation in this class. Therefore, it obviously satisfies 
all restrictions which hold for admissible transformations |15j . For example, it 
is known for a long time that for any point (and even contact) transformation 
connecting a pair of (1 + l)-dimensional evolution equations its component cor- 
responding to t depends only on t, cf. [T7J. The equations in the pair can also 
coincide. As a result, the same restriction should be satisfied by any point or con- 
tact symmetry transformation of every (1 + l)-dimensional evolution equation. 
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The simplest description of admissible transformations is obtained for normal- 
ized classes of differential equations. Roughly speaking, a class of (systems of) 
differential equations is called normalized if any admissible transformation in this 
class is induced by a transformation from its equivalence group. Different kinds of 
normalization can be defined depending on what kind of equivalence group (point, 
contact, usual, generalized, extended, etc.) is considered. Thus, the usual equiv- 
alence group G~ of the class C\s consists of those point transformations in the 
space of variables and arbitrary elements, which are projectable on the variable 
space and preserve the whole class C\s- The class C\s is called normalized in the 
usual sense if the set T(C\s) is generated by the usual equivalence group G~. As a 
consequence, all generalizations of the equivalence group within the framework of 
point transformations are trivial for this class. See [25] for precise definitions and 
further explanations. If the class C\s is normalized in certain sense with respect 
to point transformations, the point symmetry group Gg of any equation Cg from 
this class is contained in the projection of the corresponding equivalence group 
of C\s to the space of independent and dependent variables (taken for the value 
6 = 9q in the case when the generalized equivalence group is considered). 

As a rule, calculations of certain common restrictions on admissible trans- 
formations of the entire normalized class or its normalized subclasses or point 
symmetry transformations of a single equation from this class have the same level 
of complexity. For example, in order to derive the restriction that the transfor- 
mation component corresponding to t depends only on t, we should carry out 
approximately the same operations, independently of considering the whole class 
of (1 + l)-dimensional evolution equations, any well-defined subclass from this 
class or any single evolution equation. This is why it is worthwhile to first con- 
struct nested series of normalized classes of differential equations by starting from 
a quite general, obviously normalized class, imposing on each step additional aux- 
iliary conditions on the arbitrary elements and then studying the complete point 
symmetries of a single equation from the narrowest class of the constructed series. 

In the way outlined above we have already investigated hierarchies of normal- 
ized classes of generalized nonlinear Schrodinger equations [25], (l + l)-dimensional 
linear evolution equations [26], (l + l)-dimensional third-order evolution equations 
including variable-coefficient Korteweg-de Vries and modified Korteweg-de Vries 
equations [27] and generalized vorticity equations arising in the study of local 
parameterization schemes for the barotropic vorticity equation |23j . 

If an equation does not belong to a class whose admissible transformations have 
been studied earlier, one can try to map this equation using a point transformation 
to an equation from a class for which constraints on its admissible transformations 
are known a priori. Then one can either map the known constraints on admissible 
transformations back and then complete the calculations of point symmetries of 
the initial equation using the direct method or calculate the point symmetry group 
of the mapped equation using the direct method and then map this group back. 
The example on the application of this trick to the barotropic vorticity equation 
in presented in Section 01 
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3 Calculations based on Lie invariance algebra 
of the barotropic vorticity equation 



The barotropic vorticity equation on the /3-plane reads 



Ct + ipxCy ~ i>y(x + filpx = 0, 



(1) 



where ip = ip(t,x,y) is the stream function and £ := ip xx + ip yy is the relative 
vorticity, which is the vertical component of the vorticity vector. The barotropic 
vorticity equation in the formulation ((TJ is valid in situations where the two- 
dimensional wind field can be regarded as almost non-divergent and the motion in 
North-South direction is confined to a relatively small region. It is then convenient 
to use a local Cartesian coordinate system. In such a coordinate system, the effect 
of the sphericity of the Earth is conveniently taken into account by approximating 
the normal component of the vorticity due to the rotation of the Earth, 2Q sin (p, 
by its linear Taylor series expansion, where Q is the angular rotation of the Earth 
and tp is the geographic latitude. This linear approximation at some reference 
latitude ipo is given by 20 sin (^o + fly, where (3 = 2Qcosip/a and a is the radius 
of the Earth. This is the traditional /3-plane approximation, see [22] for further 
details. Then, taking the vertical component of the curl of the two-dimensional 
ideal Euler equations and using the /3-plane approximation leads to Eq. ([I]). 

It is straightforward to determine the maximal Lie invariance algebra Q of 
Eq. ([!]) using infinitesimal techniques: 



g=(V,d t ,d y ,X(f),Z(g)), 
where V = td t - xd x - yd y - 3tpd^, X(f) = f(t)d x - ft{t)yd^ and Z(g) = g(t)d^, 



and / and g run through the space of smooth functions of t. (In fact, the precise 
interpretation of q as a Lie algebra strongly depends on what space of smooth 
functions is chosen for / and g, cf. Note A.l in [8l p. 178].) This result was 
first obtained in [13] and is now easily accessible in the handbook [121 P- 223]. 
See also [2] for related discussions and the exhaustive study of the classical Lie 
reductions of Eq. ([1]). 

The nonzero commutation relations of the algebra q in the above basis are 
exhausted by the following ones: 



It is easy to see from the commutation relations that the Lie algebra q is solvable 
since 



[dt,V]=dt, [dy,V] = -dy, 

[V, X(f)\ = X(tf t + f), [V, 2(g)] = Z(tg t + 3g), 
[d t ,X(f)]=X(f t ), [d t ,Z(g)]=Z(g t ), [d y ,X(f)} 



Z(ft)- 



0' = [Q,B] = (dt,d y ,X(f),Z(g)) 
g" = [ Q ',g'} = (X(f),Z(g)), 
W = [0",0"]=O. 
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Therefore, the radical r of g coincides with the entire algebra g. The nil-radical 
of g is the ideal 

n=(d y ,X(f),Z(g)). 

Indeed, this ideal is a nilpotent subalgebra of q since 

n ( 2 )= n '=[n,n] = (Z( 5 )), n( 3 ) = [n, n'] = 0. 

It can be extended to a larger ideal of q only with two sets of elements, {dt} and 
{V, dt}- Both resulting ideals are not nilpotent. In other words, n is the maximal 
nilpotent ideal. 

Continuous point symmetries of Eq. ([1]) are determined from the elements of q 
by integration of the associated Cauchy problems. It is obvious that Eq. ([T]) also 
possesses two discrete symmetries, (t,x,y,i(j) i-> (—t,—x,y,ip) and (t,x,y,tp) t-t 
(t,x, —y, —ip), which are independent up to their composition and their composi- 
tions with continuous symmetries. The proof that the above symmetries generate 
the entire point symmetry group was, however, outstanding. 

Theorem 1. The complete point symmetry group of the barotropic vorticity equa- 
tion on the j3 -plane (pQ) is formed by the transformations 

T: t = T 1 t + T , x=—x + f(t), y = ^y + Y , 

where Ti ^ 0, e = ±1 and f and g are arbitrary functions oft. 

Proof. The discrete symmetries of the barotropic vorticity equation on the (3- 
plane are computed as described in section 12.11 The general form of a point 
transformation of the vorticity equation is: 

T: (t,x,y,4>) = (T,X,Y,V), 

where T, X, Y and ^ are regarded as functions of i, x, y and ip, whose joint 
Jacobian does not vanish. To obtain the constrained form of T, we use the above 
four proper nested megaideals of q, namely n', q" , n and q' , and q itself. Recall once 
more that the transformation T must satisfy the conditions %n' = n', %g" = g", 
%n = n, T*q' = q' and %Q = g in order to qualify as a point symmetry of the 
vorticity equation, where % denotes the push-forward of T to vector fields. In 
other words, we have 

%Z(g) = g{T^d~ t + X^d- X + Y^dy + vi/^) = Zip), (2) 
%X{f) = X{ff) + Z(gf), (3) 
%d t = T t 8- t + X t dr x + Y t 8y + Vtdj = ai d- t + a 2 d y + X{f) + Z{g\ (4) 

T.dy = Tydl + Xyd X + YyOy + Vydj = b 1 8y + X(f) + Z^) , (5) 

%V = C1 V + c 2 d~ t + c 3 8y + X(f D ) + Z(g D ), (6) 
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where all /'s and g's are smooth functions of t which are determined, as the 
constant parameters a%, 0,2, fei, c±, 02 and C3, by % and the operator from the 
corresponding left-hand side. 

We will derive constraints on 7*, consequently equating coefficients of vector 
fields in conditions ©-([B]) and taking into account constraints obtained on previ- 
ous steps. Thus Eq. ([2]) immediately implies = = = (hence ^ / 0) 
and g^^ = g 9 . Evaluating the last equation for g = 1 and g = t and combin- 
ing the results gives t = (T)/^ 1 (T) , where g l = g 9 \ g =t arid g 1 = g 9 \ g =i- As the 
derivative with respect to T in the right hand side of this equality does not vanish, 
the condition T = T(t) must hold. This implies that depends only on t. 

As then %X(f) = fX x 8 x + fY x 8 y + (f* x - ftyV^dj, it follows from Eq. © 
that Y x = and 

fX x = p , fV x - f t y^ =-ffY + gf. 

Evaluating the first of the displayed equalities for / = 1, we derive that X x = 
/ 1 (T) =: Xi(t). Therefore, /* (T) = f(t)Xi(t). The second equality then reads 

f^-fty^ = - ( -^Y + gf. 

J-t 

Setting / = 1 and / = t in the last equality and combining the resulting equalities 
yields j/*^ = (T t )~ 1 XiY +tg x -~g l , where g l = ~gf\ f=t and g 1 = gf\ f=1 . As X l ^ 
this equation implies that Y = Y\{t)y + Yo(t). 

After analyzing Eq. (j4]), we find Tt = const, Yf = const, which leads to 
Y\ = const, X t = f(T) and thus X tx = 0, i.e., X\ = const. Finally, Eq. Q 
also implies = —fjY + 5- In a similar manner, upon taking into account 
the restrictions already derived so far, collecting coefficients in Eq. ((5]) gives the 
constraint X y = f v =: X2 = const since X yt = 0. Moreover, ^ y = g y , as f~ = 0. 

The final restrictions on T based on the preservation of 3 are derivable from 
Eq. dl]), where 

%V = tT t d- t + {tX t - xX x - yX y )3 x + (tY t - yY y )d y 
+ (i¥ t - xV x - y^! y - 3^)^- 

Collecting the coefficients of <9j and d y , we obtain that c\ = 1 and Y t = 0. Simi- 
larly, equating the coefficients of dy and further splitting with respect to x implies 
that ty x = 0. 

The results obtained so far lead to the following constrained form of the general 
point symmetry transformation of the vorticity equation ([TJ) 

T = T 1 t + T , X = X 1 x + X 2 y + f(t), Y = Y iy + Y , 
* = * 1 ^ + * 2 (*)y + *4(*), 

where To, T\, X\, X2, Yq, Y\ and ^1 are arbitrary constants, T\X{Yi^i ^ 0, 
and f(t), *&2(t) an d ^(i) are arbitrary time-dependent functions. The form ([7]) 
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takes into account all constraints on point symmetries of ([I]), which follow from 
the preservation of the maximal Lie invariance algebra q by the associated push- 
forward of vector fields. 

Now the direct method should be applied. We carry out a transformation of the 
form ((7J) in the vorticity equation. For this aim, we calculate the transformation 
rules for the partial derivative operators: 

1 (a ft a \ a 1 Q Q 1 ( ^ X 2 



Further restrictions on T can be imposed upon noting that the term iptxy can 
only arise in the expression for -0^--- , which is 

2*i X 2 

^tyy ~ ~TiY 1 X 1 V ^ ' 

This obviously implies that X 2 = 0. In a similar fashion, the expression for Q is 

^ = f\ ( (Xtf Ct + (jYtf ~ (X{f) ijyyt ) + " ' ' 

upon using ip xx t = Ct — ^yyt- Hence (^i) 2 = (Y\) 2 as there are no other terms 
with 'ipyyt in the invariance condition. After taking into account these two more 
restrictions on T, it is straightforward to expand the transformed version of the 
vorticity equation. This yields 

*1 r /■ , (*l) 2 , /- f*l , , *5 



T 1 (X 1 )^ L T X {X 
'—ib — 



! , (*l) 2 , *2 \ *1 , 

l)3 ^+ {x{fYS xt * v \Y X % + Yj (X^ 



+ fi—^x = ^ fvA2 (Ct + ^xCy ~ ^yCx + P^x 



The invariance condition is fulfilled provided that the constraints 

hold. This completes the proof of the theorem. ■ 

Corollary 3. The barotropic vorticity equation on the /3-plane possesses only two 
independent discrete point symmetries, which are given by 

Ti : (t, x, y, ip) ^ (-*, -x, y, ip), F 2 : (t, x, y, ip) i-> (i, x, -y, -ip). 

They generate the group of discrete symmetry transformations of the barotropic 
vorticity equation on the (3-plane, which is isomorphic to Z 2 xZ 2 , where 1? denotes 
the cyclic group of two elements. 
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4 Direct method and admissible transformations 
of classes of generalized vorticity equations 

The construction of the complete point symmetry group G of the barotropic vor- 
ticity equation ([1]) by means of using only the direct method involves cumbersome 
and sophisticated calculations. As Eq. ([I]) is a third-order PDE in three indepen- 
dent variables, the system of determining equations for transformations from G 
is an overdetermined nonlinear system of PDEs in four independent variables, 
which should be solved by taking into account the nonsingularity condition of 
the point transformations. This is an extremely challenging task. Fortunately, a 
hierarchy of normalized classes of generalized vorticity equations was recently con- 
structed |23j that allows us to strongly simplify the whole investigation. Eq. ([1]) 
belongs only to the narrowest class of this hierarchy, which is quite wide and 
consists of equations of the general form 

Ct = F(t, X, y, 1p, Ipx-ilpyi Cj Cxi Cy: Cxxi Cxyi Cyy)i C := ^xx ~\~ ^yyt (8) 

where (F Cx , F Cy , F Cxx , F (xy , F (yy ) / (0, 0, 0, 0, 0) . The equivalence group G? of this 
class is formed by the transformations 

t = T{t), x = Z 1 {t,x,y), y = Z 2 (t,x,y), $ = T(t)^ + x, y), 

^4(r-(rV + (x) -?(rO + ©/ + m/ 

where T, Z 1 , T and $ are arbitrary smooth functions of their arguments, satisfying 
the conditions Z\Z\ = 0, Z\Z\ = Z\Z\ := L and T t TL ^ 0. The subscripts 1 
and 2 denote differentiation with respect to x and y, respectively, the indices i 
and j run through the set {1,2} and the summation over repeated indices is 
understood. As Eq. ([T]) is an element of the class ([8]) and this class is normalized, 
the point symmetry group G of Eq. ([1]) is contained in the projection G± of the 
equivalence group G± of the class ([8]) to the variable space (t,x,y,ip). At the 
same time, the group G is much narrower than the group Gi, and in order to 
single out G from G± we should still derive and solve a quite cumbersome system 
of additional constraints. Instead of this we use the trick described in the end of 
Section 12.21 Namely, by the transformation 

V> = ^ + ^ 3 , (9) 

which identically acts on the independent variables and which is prolonged to the 
vorticity according to the formula ( = C + f3y, we map Eq. ([T]) to the equation 

R 

Ct + ^xly ~ IpyQx = -^V 2 Cx- (10) 
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Eq. (jlOj) belongs to the subclass of class ([8]) that is singled out by the constraints 
= 0, = 0, F^ x = —Q y and = Cx, he., the class consisting of the equations 
of the form 

Ct ~\~ 'tpxCy ~ IpyCx = H(t, X, y, (^ x , Qy, Qxxi Cxyi Cyy)i C := ^xx ~\~ ^yyi (H) 

where H is an arbitrary smooth function of its arguments, which is assumed as 
an arbitrary element instead of F = H — if) x Cy + VfyCz- The class fill j) also is a 
member of the above hierarchy of normalized classes. Its equivalence group G% 
is much narrower than Gi and is formed by the transformations 

1 ~ w , „.„\ i 2 



t = r, x = X(xc — ys) + 7 , ey = \(xs + yc) + 7 
rp = e— [ \ib 

n 

e ( „ A+ . > , „ \ S v +ay _ 5 x +(?x _ 2 / a 



+ ^0 t (x 2 +y 2 ) - jj(xs+yc) + 7 2 (xc-ys)^ + 5 + |(x 2 +y 2 ), 



where e = ±1, c = cos 6,5 = s'mO; r, A, 9, 7* and a are arbitrary smooth functions 
of t satisfying the conditions A > 0, t u = and r t ^ and 5 = 5(t,x,y) runs 
through the set of solutions of the Laplace equation S xx + 5 yy = 0. 

In order to derive the additional constraints that are satisfied by the group 
parameters of transformations from the point symmetry group G2 of Eq. (|10p , we 
substitute the values H = —(3y 2 (,x/2 and H = —f3y 2 ( iX /2 as well as expressions 
for the transformed variables and derivatives via the initial ones into the trans- 
formation component for H and then make all possible splitting in the obtained 
equality. As a result, we derive the additional constraints 



2 



1 e/? 7 2 , , , £/?(7 2 ) 2 

-, ^ = ^TT> S x = -ax, 5 y = ay + ^-. 



After projecting transformations from on the variable space (t,x,y,ip), con- 
straining the group parameters using the above conditions and taking the adjoint 
action of the inverse of the transformation (J9|), we obtain, up to re-denoting, the 
transformations from Theorem [TJ 



5 Conclusion 

In this paper, we have computed the complete point symmetry group of the 
barotropic vorticity equation on the /3-plane. It is obvious that both of the tech- 
niques presented in this paper are applicable to general systems of differential 
equations. 

Despite of the apparent simplicity of the techniques employed above, there are 
a number of features that should be discussed properly. In particular, the relation 
between discrete symmetries of a differential equation and discrete automorphisms 
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of the corresponding maximal Lie invar iance algebra is neither injective nor sur- 
jective. This is why it can be misleading to restrict the consideration to discrete 
automorphism when trying to finding discrete symmetries. This and related issues 
will be investigated and discussed more thoroughly in a forthcoming work. 
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